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The Ising model in clustered scale-free networks has been studied by Monte Carlo simulations. 
These networks are characterized by a degree distribution of the form P{k) ~ k ^ for large k. 
Clustering is introduced in the networks by inserting triangles, i.e., triads of connected nodes. The 
transition from a ferromagnetic (FM) to a paramagnetic (PM) phase has been studied as a function 
of the exponent 7 and the triangle density. For 7 > 3 our results are in line with earlier simulations, 
and a phase transition appears at a temperature Tc( 7 ) in the thermodynamic limit (system size 
N —>■ 00 ). For 7 < 3, a FM-PM crossover appears at a size-dependent temperature Tco, so that the 
system remains in a FM state at any finite temperature in the limit N 00 . Thus, for 7 = 3, Too 
scales as InA^, whereas for 7 < 3, we find Too ~ JN^, where the exponent 2 decreases for increasing 
7 . Adding motifs (triangles in our case) to the networks causes an increase in the transition (or 
crossover) temperature for exponent 7 > 3 (or < 3). For 7 > 3, this increase is due to changes in 
the mean values {k) and (fc^), i.e., the transition is controlled by the degree distribution (nearest 
neighbor connectivities). For 7 < 3, however, we find that clustered and unclustered networks 
with the same size and distribution P{k) have different crossover temperature, i.e., clustering favors 
FM correlations, thus increasing the temperature Tco. The effect of a degree cutoff fccut on the 
asymptotic behavior of Tco is discussed. 


I. INTRODUCTION 

Many natural and artificial systems have a network 
structure, with nodes representing typical system units 
and edges playing the role of interactions between con¬ 
nected pairs of units. Complex networks can be used to 
model various kinds of real-life systems (social, economic, 
technological, biolc^ical), and to analyze processes tak¬ 
ing place on them [iH^- In recent years, various network 
models have been designed capturing aspects of real sys¬ 
tems, thus allowing to explain empirical data in several 
fields. This is the case of small-world @ and scale-free 
networks Q , which provide us with the underlying topo¬ 
logical structure to analyze processes such as signal prop- 
agation [6|,H , as well as the spread of information [^, , 

opinions [lll| and infections [l^ . These types of net¬ 
works have been also used to study statistical physical 
problems as percolation [ll,[i3| and cooperative phenom¬ 
ena [iB-[2i|. 

In scale-free (SF) networks the degree distribution 
P(fc), where k is the number of links connected to a node, 
has a power-law decay Psf{k) ^ k~'^ [H, This type 
of networks have been found in several real-life systems, 
such as the internet [2^, the world-wide web [l^, pro¬ 
tein interaction networks [2^, and social systems [2^ . 
In both natural and artificial systems, the exponent 7 
controlling the deg ree distribution is usually in the range 
2 < 7 < 3 m, HH. The origin of power-law deme dis¬ 
tributions was studied by Barabasi and Albert Q, who 
found that two ingredients can explain the scale-free na¬ 
ture of networks, namely growth and preferential attach¬ 
ment. More general models based on these ingredients 
have appeared later in the literature [H, [2^ . One can 
also deal with equilibrium SF networks, defined as statis¬ 
tical ensembles of random networks with a given degree 


distribution Psf(fc), for which one may analyze several 
properties as a function of the exponent 7 [^, . 

Many real-life networks include clustering, i.e., the 
probability of finding loops of small size is larger than 
in random networks. This has been in particular quan¬ 
tified by the so-called clustering coefficient, which mea¬ 
sures the likelihood of “triangles” in a network Q . Most 
network models employed in the past did not include 
clustering. Some of them, such as the Watts-Strogatz 
small-world model @ show clustering, but are not well- 
suited as models of most actual networks. Several compu¬ 
tational models of clustered networks have been defined 
along the years [sil - ls^ . but in general their properties 
cannot be calculated by analytical procedures. In last 
years, it was shown that generalized random graphs can 
be generated incorporating clustering in such a way that 
exact formulas can be derived for many of their proper¬ 
ties. This is the case of the networks defined by Newman 
[13 and Miller [H. 

Cooperative phenomena in complex networks are 
known to display characteristics related to the partic¬ 
ular topology of these systems M- The Ising model on 
SF networks has been studied by using several theoretical 
techniques and its critical behavior was found 

to depend on the exponent 7. Two different regimes 
appear for uncorrelated networks. On one side, for an 
exponent 7 > 3, the average value {k'^) is finite in the 
large-size limit, and there appears a ferromagnetic (FM) 
to paramagnetic (PM) transition at a finite temperature 
Tc- On the other side, when (fc^) diverges (as happens 
for 7 < 3), the system remains in its ordered FM phase 
at any temperature, so that there is no phase transition 
in the thermodynamic limit. The antiferromagnetic Ising 
model has been also studied in scale-free networks, where 
spin-glass phases have been found [4l|, [13 ■ 
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All this refers to unclustered random networks with 
a power-law degree distribution. One may ask how this 
picture changes when the networks are clustered, i.e., the 
clustering coefficient has a non-negligible value. In prin¬ 
ciple, one expects that the presence of small loops in the 
networks will enhance correlations between spins located 
on network nodes, thus favoring ordered schemes such as 
an FM pattern. Thus, the effects of clustering on various 
cooperative phenomena in complex networks have been 
studied earlier. This is the case of percolation 1^ 
epidemics [M 113) mi ) and dynamical processes [is 
Yoon et al. ^ studied the Ising model on networks with 
arbitrary distribution of motifs using an analytical proce¬ 
dure, the so-called belief-propagation algorithm. For the 
networks considered in that paper, where the thermody¬ 
namic limit is well defined (Tc converges for N —>■ oo), 
these authors found that clustering increases the critical 
temperature in comparison with tree-like networks with 
the same mean degree, but does not change the critical 
behavior. 

In this paper we study the FM-PM transition for the 
Ising model in scale-free networks with clustering, which 
is realized by introducing triangles in the networks, i.e., 
three-membered loops. Several values of the exponent 7 
are considered, as well as various concentrations of tri¬ 
angles. We employ Monte Carlo (MC) simulations to 
obtain the transition temperature, when it is well de¬ 
fined (for 7 > 3 ), and to derive the size dependence of 
the crossover temperature for cases where (A:^) diverges 
as Y —>■ 00 (7 < 3). 

The paper is organized as follows. In Sec. II we de¬ 
scribe the clustered networks considered here. In Sec. Ill 
we present the computational method employed to carry 
out MC simulations. In Sec. IV we present results of the 
simulations and a discussion for the different parameter 
regions, depending on the value of the exponent 7 (> 3, 
= 3, or < 3). The paper closes with a summary in Sec. V. 


II. SCALE-FREE NETWORKS WITH 
CLUSTERING 

We consider clustered networks with a degree dis¬ 
tribution P{k) that follows the power-law dependence 
Psi{k) ^ k~^ for large degree k. Clustering is included 
by inserting triangles in the networks, i.e., triads of con¬ 
nected nodes. Other kinds of polygons (squares, pen¬ 
tagons, ...) can be introduced to study their effect on 
critical phenomena in physical systems, but we choose 
triangles since they cause stronger correlations between 
entities defined on network sites, as in the case of the 
Ising model considered here. 

We generate networks by following the method pro¬ 
posed by Newman Q, where one separately specifies 
the number of edges and the number of triangles. This 
procedure allows one to generalize random graphs to in¬ 
corporate clustering in a simple way, so that exact for¬ 
mulas can be derived for many properties of the resulting 



FIG. 1: Schematic representation of a typical network con¬ 
sidered in this work, for which one separately specifies the 
number of single links and triangles attached to each node. 
Triangles are indicated by asterisks (*). 


networks [s^. 

For a network of size N (number of nodes), we call 
ti the number of triangles in which node i takes part, 
and Si the number of single edges not included in the 
triangles. This means that edges within the triangles 
are listed separately from single links. Thus, a single 
link can be viewed as a network element joining together 
two nodes and a triangle as an element connecting three 
nodes. The degree ki of node i is then ki = st + 2ti, as 
each triangle connects it to two other nodes. A picture 
of such a network is presented in Fig. 1, where triangles 
are indicated with asterisks (*). 

To generate the networks, we first define the edges. 
We assign to each node i a random number s^, which 
represents the number of outgoing links from this node 
(stubs). The set of numbers {si}A ^ (with Si > ko, the 
minimum allowed degree) is taken from the probability 
distribution Psf{s) ^ [s^, giving a total number of 

stubs K = Xlfci We impose the restriction that K 
must be an even integer. Then, we connect stubs at 
random (giving a total of L = Ar /2 connections), with 
the conditions: (i) no two nodes can have more than one 
bond connecting them (no multiple connections), and (ii) 
no node can be connected by a link to itself (no self¬ 
connections). 

In a second step we introduce triangles into the net¬ 
works. Their number Ya is controlled by the parameter 
v, which gives the mean number of triangles in which 
a generic node is included (Ya = The number 

of triangles ti associated to a node i is drawn from a 
Poisson distribution Q{t) = jt]. Thus, we have ti 

“corners” associated to node i, and the total number is 
T = = 3Ya. We impose the condition that 

T be a multiple of 3. Then, we take triads of corners 
uniformly at random to form triangles, taking into ac¬ 
count conditions (i) and (ii) above to avoid multiple and 
self-connections. Note that single links can by chance 
form triangles. Calling Nt the number of such triangles. 
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FIG. 2: Probability density as a function of the degree k 
for networks with 7 = 3 , minimum degree ko = 3, and N = 
10® nodes. The data shown are an average over 200 network 
realizations for each value of the parameter v — 0, 1, and 2. 


their density Nt/N for a given mean degree (k) vanishes 
as Ai —>■ 00 . In fact, Nt/N scales as 1/N for large N 

SUllIil- 

In complex networks, one usually defines the clustering 
coefficient C as the ratio C = SNa/N^, where N 3 is the 
number of connected triplets [^: 

N^ = NY,^^^Y^P{k) = lN{{k^)-{k)) . ( 1 ) 

k 

Thus, for the networks discussed here, we have 

and the clustering coefficient can be changed as a function 
of the parameter v. 

Aside from 7 and the triangle density v, our networks 
are defined by the minimum degree kg. Since we are in¬ 
terested in finite-size effects, the network size N is also an 
important variable in our discussion. The degree distri¬ 
bution P{k) obtained for networks generated by following 
the procedure described above is presented in Fig. 2. In 
this figure, we have plotted P{k) for networks including 
10® nodes, with 7 = 3 and ko = 3. Each curve corre¬ 
sponds to a particular value of the parameter 1 / = 0, 1, 
and 2 , including in each case an average over 200 net¬ 
work realizations. Comparing the curves corresponding 
to different ly values, one observes that the introduction 
of triangles in the networks causes clear changes in the 
distribution P{k) for small k values. However, for large 
degrees, the distribution is found to follow the depen¬ 
dence Psf(fc) ~ k~'^ typical of scale-free networks, with 
7 = 3 in the present case. This could be expected, since 
the Poisson distribution Q{t) associated to the triangles 



FIG. 3: Mean degree {k) as a function of for networks 

with 7 = 3, fco = 3, and three values of iz: 0 (squares), 1 
(circles), and 2 (diamonds). Error bars of the simulations 
results are less than the symbol size. Lines are fits to the 
expression (k) = a -|- hN~^P. 

has a much faster exponential-like decay for large t. In 
the three cases shown in Fig. 2 there appears an effective 
cutoff kcut ^ 300, which is related to the finite size N 
of the networks (see below). We note that a maximum 
degree fccut was explicitly introduced earlier in scale-free 
networks for computational convenience [s^. 

An important characteristic of the considered net¬ 
works, that will be employed below to discuss the results 
of the Ising model, is the mean degree {k). For scale-free 
networks with ^ = 0, the mean degree is given by 

00 ^ 

{k)^ = y2kP,f{k)^ko^^, (3) 

k=ko ' 

where the last expression is obtained by replacing the 
sum by an integral, which is justified for large N. Note 
that we assume here 7 > 2 and that the distribution 
Psf(fc) is normalized to unity (for 7 < 2 the mean degree 
(k) diverges in the large-size limit: {k)ao 00 ). Then, 
for our networks including triangles, one has 

(k)^ Ki 2iy + ko- —i. (4) 

7-2 

For finite networks, a size effect is expected to appear 
in the mean degree, as a consequence of the effective cut¬ 
off appearing in the degree distribution (see Fig. 2). This 
is shown in Fig. 3, where we present (k) vs for our 

generated networks with 7 = 3 and fco = 3. The mean 
degree decreases as N~^/^ increases, i.e., (k) increases as 
the system size is raised. The data shown in Fig. 3 for 
iz = 0, 1, and 2, follow a linear dependence, and in fact 
can be fitted as (fc) = a + bN~^^‘^, at least in the region 
plotted in the figure {N > 500). The fit parameter a is 
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close to the mean degree (fc)oo given in Eq. (jd)), and the 
small difference is mainly due to the replacement of sums 
by integrals in the derivation of that equation. Note that 
for 7 = 3, = 0 and fco = 3, Eq. (|1]) yields (fc)oo = 

6 . The slope b obtained from the linear fits turns out to 
be the same (within statistical noise) for the three cases 
shown in Fig. 3. 

This size dependence of (fc) can be understood by not¬ 
ing that the effective cutoff k^ut appearing in a power-law 
degree distribution is related with the network size N by 
the expression [^, 

OO 

= (5) 

^cut 


where c is a constant on the order of unity. From this 
expression, one can derive for 7 = 3 (see the Appendix): 


(k) 



Comparing with the fit shown in Fig. 3, one finds c ~ 
7, which introduced into Eq. (IAtI) yields for the cutoff 
^cut ~ 360, in line with the results shown in Fig. 2, thus 
providing a consistency check for our arguments. 

For scale-free networks with an exponent 7 < 3, Catan- 
zaro et al. found that appreciable correlations ap¬ 
pear between degrees of adjacent nodes when no multi¬ 
ple and self-connections are allowed. Such degree corre¬ 
lations can be avoided by assuming a cutoff fccut 
Thus, for 7 < 3 we generate here networks with cutoff 
fccut = For the clustered networks considered here, 

generated as in Ref. [s^, the presence of triangles intro¬ 
duces degree correlations between nodes forming part of 
a triangle in the network. 


III. SIMULATION METHOD 

On the networks defined in Sec. II, we consider a spin 
model given by the Hamiltonian 

H = -Y,.h,S,S,, (7) 

i<j 

where Si = ±1 (i = 1,..., N) are Ising spin variables, and 
the coupling matrix Jy is given by 


N. This allows us to investigate the transition from a FM 
((M) ^ 0) to a PM ((M) = 0) regime as T is increased. 

Depending on the value of the exponent 7 defining the 
power-law distribution of single edges, two different cases 
are found [H, [H, First, for 7 > 3, one expects a 
phase transition with a well-defined transition tempera¬ 
ture Tc (< 00 ) in the thermodynamic limit N —>■ 00 . Sec¬ 
ond, for 7 < 3 a FM-PM crossover is known to appear for 
scale-free networks, with a crossover temperature Tco{N) 
increasing with system size and diverging to infinity as 
N —>■ OO. 

For the cases where a FM-PM transition occurs in the 
thermodynamic limit (7 > 3), the transition temperature 
Tc has been obtained here by using Binder’s fourth-order 
cumulant 


Un{T) = I - 


{M^)n 


( 9 ) 


The average values in this expression are taken over dif¬ 
ferent network realizations and over different spin con¬ 
figurations for a given network at temperature T. In 
this case, the transition temperature is obtained from the 
unique crossing point of the functions Un{T) for several 
system sizes N j^, . 

In the second case (7 < 3), the size-dependent 
crossover temperature Tco{N) has been obtained from 
the maximum of the magnetization fluctuations (AM)^ 
as a function of temperature, with 


iAM)% = {M^)n - {M)% . (10) 


We note that Tco values derived by using this criterion 
agree within error bars with those found from the maxi¬ 
mum derivative of the heat capacity . 

The largest networks considered here included about 
10® sites. Such network sizes were employed in particu¬ 
lar to study the dependence of Tea on N for 7 < 3. For 
the cases where a phase transition exists in the thermo¬ 
dynamic limit (7 > 3), sizes around 4 x 10^ nodes were 
considered. The results presented below were obtained 
by averaging in each case over 800 networks, except for 
the largest system sizes, for which 400 network realiza¬ 
tions were generated. Similar MC simulations have been 
carried out earlier to study ferromagnetic [13, and 
antiferromagnetic [13, Ising models in complex net¬ 
works. 


IV. RESULTS AND DISCUSSION 


f J(>0), if f and j are connected, Case 7 > 3 

1 0 , otherwise. ' 


This model has been studied by means of Monte Carlo 
simulations, sampling the configuration space by using 
the Metropolis local update algorithm [sj. We are par¬ 
ticularly interested in the behavior of the magnetization 
M = '^f^iSi/N. For a given set of parameters ( 7 , z/, 
fcp) defining the networks, the average value {M) has been 
studied as a function of temperature T and system size 


For unclustered scale-free networks with an exponent 
7 > 3, the average value (fc^) converges to a finite value 
as N ^ OO. In this case, anal ytic al calculations [H, 
and Monte Carlo simulations predict a well-defined 
FM-PM transition temperature Tc given by 


J 

¥c 



(P) 


(fc2) - 2(fc) 


( 11 ) 
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FIG. 4: Transition temperature for networks with 7 = 
5 as a function of the minimum degree ko. Open symbols 
represent results of MC simulations, as obtained from Binder’s 
cumulant for v = Q (squares), 1 (circles), and 2 (diamonds). 
Error bars are less than the symbol size. Solid lines were 
derived from the analytical expression given in Eq. (HU, using 
values of (fc) and (fe^) obtained from Eqs. ([4]) and (I14II . Solid 
squares are data obtained in Ref. for unclustered networks 
{u = Q). The dashed line is a guide to the eye. 

This equation is equivalent to Eq. (56) in Ref. m 

For the clustered networks considered here with 7 > 3, 
we have calculated from the Binder’s cumulant Un 
for several values of the triangle density v and minimum 
degree fcp > 1. In each case, four different network sizes 
were considered. In all these cases, a well defined tran¬ 
sition temperature was found from the crossing poin t of 
the curves Un{T) for different sizes Ai, as in Ref. We 
note that (A:^) — 2 (k) > 0 for /cq > 1 , and Tc is well de¬ 
fined by Eq. m ■ For fco = 1, the simulated networks 
consist of many disconnected components, and Binder’s 
cumulant Un{T) does not give a unique crossing point 
for different system sizes N. 

Going to the results of the present MC simulations, 
in Fig. 4 we present the transition temperature Tc as 
a function of the minimum degree fco for an exponent 
7 = 5, and three values of the triangle density v = 0, 
1, and 2. In the three cases we observe a linear increase 
of Tc for rising fcp. The case 1 ^ = 0 corresponds to a 
power-law degree distribution Psf(fc) ~ fc~^ (unclustered 
networks). Values of Tc found here for these networks are 
somewhat higher than those obtained in Ref. by an 
amount ATc ^ 0.6 J, due to the strict degree cutoff fccut 
employed in that work. These earlier results are shown 
in Fig. 4 as open symbols. For 1 / = 1 we find Tc values 
higher than for 1 / = 0 , and the transition temperature 
increases further, by the same amount, for v = 2. 

Looking at Eq. m, one expects that Tc should have 
an explicit dependence on (fc) and (fc^), rather than the 
minimum degree fco itself. In Fig. 5 we present the tran¬ 



FIG. 5: Transition temperature for scale-free networks 
with 7 = 5 as a function of the mean degree (fc). Data points 
are results for networks with = 0 (open squares) and v 
= 2 (solid diamonds). Error bars are less than the symbol 
size. Solid and dashed-dotted lines are analytical predictions 
from Eq. m for iz = 0 and 2, respectively. The dashed line 
represents the mean-field result given by Eq. m- 

sition temperature Tc as a function of the mean degree 
(fc), as derived from MC simulations for unclustered net¬ 
works (y = 0 , open squares) and clustered networks with 
V = 2 (solid diamonds). Lines in this figure were ob¬ 
tained from the average values (fc) and (fc^). Taking into 
account that ki = Si + 2ti, the average value (fc^) for 
clustered networks with 7 > 3 can be calculated as: 

(fc 2 ) = (s2)+4(s)(f)+4(t2), ( 12 ) 

since s and t are independent due to the way of building 
up these networks. For the power-law distribution of Si 
corresponding to single edges, we have in the large-V 
limit: 

(s^)oo~fco^^—^ (13) 

7-3 

so that 

(fc^)oo = fco-7 + 4 fco iz- - -|- 4 ^ (jz -|- 1) . (14) 

7—3 7“2 

Here v and v {v + 1) are the average values (t) and (t^) 
corresponding to the Poisson distribution Q{t) of trian¬ 
gles in these networks. Introducing Eq. (11^1) for (fc^)oo 
and Eq. dH) for (fc)oo into Eq. (fTTI) . with the parameters 
7 = 5 and iz = 0, we find for the transition temperature 
Tc the solid line shown in Fig. 5. This line lies very close 
to the results derived from our Monte Carlo simulations 
for unclustered networks (open squares). Similarly, for 
7 = 5 and ly = 2, we obtain from Eq. (HID the dashed- 
dotted line, which coincides with the data obtained from 
simulations for clustered networks (full diamonds). Note 
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that for 7 > 3, finite-size effects on (k) and (fc^) are 
negligible for the network sizes employed in our simula¬ 
tions, so that the agreement between Eq. (HU and our 
simulation data is good. In fact, the MC results agree 
within error bars with the transition temperature given 
by Eq. (fTT|) . 

For comparison, we also present in Fig. 5 the critical 
temperature obtained in a mean-field approach [iU| : 

( 15 ) 

which is displayed as a dashed line [v = 0). Note that 
this mean-field expression can be derived from Eq. (HU 
in the limit {k^)/{k) ^ 1. Expanding Eq. (fTTl) for small 
{k)/{k'^), one has 


Tc jk^) Uk) . {k? \ 

j {k) 3{k^) \{k^)y 


(16) 


where we recognize the first term in the expansion as the 
mean-field approximation in Eq. HU- 

The critical temperature Tc derived from our MC sim¬ 
ulations for 7 = 5, and shown in Fig. 4 as a function 
of the minimum degree ko, can be fitted linearly with 
good precision as Tc = a ko + b. In fact, for = 0 we 
find a = 1.52 and b = —1.17. The value of a can be es¬ 
timated from the mean-field approximation in Eq. HU, 
which yields ~ Sfcg J/2. 

Turning to the results found for clustered networks, we 
observe in Fig. 4 that, for a given kg, the transition tem¬ 
perature clearly increases when the triangle density rises. 
However, the same expression for Tc given in Eq. HU 
reproduces well the MC results for clustered and unclus¬ 
tered networks, once the corresponding values for (fc) and 
(fc^) are introduced, as shown in Fig. 5. Since these aver¬ 
age values depend only on the degree distribution P{k) 
(nearest neighbors), this means that the transition tem¬ 
perature for networks with 7 > 3 does not depend on the 
clustering. Thus, including triangles in these networks 
changes the transition temperature because it changes 
the degree distribution P{k), but networks with the same 
P{k) but without triangles give the same Tc, as predicted 
byEq. (HU. This does not happen for the crossover tem¬ 
perature Tea obtained for finite-size networks with 7 < 3 
(see below). 


B. Case 7 = 3 


For unclustered scale-free networks with an exponent 
7 close to, but higher than 3, the transition temperature 
Tc is given from Eq. (HU by 


^ ^ r. 7 - 2 

J (fc) ° 7 — 3 


(17) 


Then, dTc/dj Ki —^ 0^/(7 — 3)^ < 0, and Tc increases 
as 7 is reduced, eventually diverging for 7 —>■ 3, as a 
consequence of the divergence of (fc^). 



FIG. 6 : Crossover temperature TcojJ for scale-free networks 
with 7 = 3 and minimum degree fco = 3, as a function of the 
system size N. Squares, v = 0; circles, v = 1; diamonds, v = 
2. Lines are least-square fits to the data points for N > 4000 
nodes. 


For 7 = 3 , analytical calculations [H, [s^ have pre¬ 
dicted a FM-PM crossover at a size-dependent tempera¬ 
ture Tco, which scales as logfV. This dependence of the 
crossover temperature agrees with that derived from MC 
simulations for the same type of networks [s^. A loga¬ 
rithmic increase of Tco with system size N has been also 
found by Aleksiejuk et al. from MC simulations 

of the Ising model in Barabasi-Albert growing networks. 
Note that these networks (with 7 = 3 ) display correla¬ 
tions between degrees of adjacent nodes [^. 

For scale-free networks with 7 = 3 and ^ = 0, the mean 
degree (fc) can be approximated as [see Eqs. dU and dU]: 


{k) 



(18) 


and (fc^) is given by [see Appendix, Eq. (IA13|) ]: 

{k^)K.kl\aN. (19) 

Applying Eq. HU to the size-dependent crossover tem¬ 
perature corresponding to 7 = 3, one finds for large sys¬ 
tem size N-. 

Tco-ifcoJlnA^ . (20) 

For uncorrelated scale-free networks with 7 = 3 , Doro- 

govtsev et al. [ 1 ^ found for the crossover temperature 

Tco-i(fc)JlnAf, (21) 

which coincides with Eq. dini) for (k) Ri fco(7 - l)/(7 - 
2 ) = 2 fco. 

In Fig. 6 we present the results of our MC simulations 
for Tco as a function of the network size N. Data points 
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correspond to networks with 7 = 8 and three values of 
the triangle density ly — 0, 1, and 2. The observed lin¬ 
ear trend of the data points in this semilogarithmic plot 
indicates a dependence Tco JhiN, as that given in 
Eq. dini). We find that values of the crossover tempera¬ 
ture for system size N < 4000 tend to be higher than the 
linear asymptotic trend found for larger sizes. In fact, in 
the linear fits presented in Fig. 6 , we only included sizes 
N > 4000. The deviation for small N is particularly 
observed for :/ = 0 and 1. Thus, our results indicate a 
dependence T^o/ J = A\n N + B , with a constant A that 
decreases for increasing v. We found A = 1.38, 1.25, and 
1.09 for V = 0, 1, and 2, respectively. Note that in the 
case v = 0 (scale-free networks without clustering), the 
slope A is somewhat smaller than the value predicted by 
Eq. (I^Dl) for ko = 3, i.e., A = ko/2 = 1.5. 

A similar logarithmic dependence of Tco upon N was 
observed in earlier works. For scale-free networks with a 
strict cutoff kcut, the prefactor A was found to increase 
linearly with kg, so that A/kg = 0.28 ± 0.01 [s^. For 
Barabasi-Albert networks with fco = 5, Aleksiejuk et al. 
[ 5 ^ found from a fit similar to ours A = 2.6, which means 
A/kg = 0.52, similar to our A/kg = 0.46 for i/ = 0 . 

Given the increase in crossover temperature for rising 
system size N, it is worthwhile analyzing the dependence 
of Tco on the minimum degree kg. In fact, given the pa¬ 
rameters 7 and v, kg controls the mean degree (k) of our 
networks. In Fig. 7(a) we display results for Tco/J for 
an exponent 7 = 3 and triangle density ly = 0, 1, and 2. 
In all cases, the networks included N = 8000 nodes. As 
expected, for a given value of v, the crossover tempera¬ 
ture increases as (fc) (or kg) is raised. Moreover, the line 
giving the dependence of Tco on (k) shifts downwards for 
rising v. In view of Eq. (HZl, this can be interpreted from 
a decrease in the ratio (fc^)/(fc) for networks with con¬ 
stant size and increasing triangle density v. Note, how¬ 
ever, that the dependence of Tco on (k) is not strictly 
linear for fixed N, as predicted by Eq. (HI}. This is a 
finite-size effect, since this equation corresponds to the 
asymptotic limit, valid in the large-A^ regime, so that for 
N = 8000 such an effect is still clearly appreciable in the 
results shown in Fig. 7(a). This can be further visualized 
in Fig. 7(b), where we present the ratio Tco/{J{k)) for 
the same data as in panel (a). Values of this ratio corre¬ 
sponding to different triangle densities v converge one to 
the other as the mean degree increases. 

For networks with a given size N, it is interesting to 
analyze the dependence of Tco on the parameter v for a 
fixed value of the mean degree (fc). Since (fc) = (s) -|- 2v, 
one can obtain networks with given (fc) by simultane¬ 
ously changing (s) and v. In the actual implementations, 
we varied the minimum degree kg, which defines a mean 
value (s), and then we took v = {{k) — (s))/2, (k) be¬ 
ing the required mean degree of the clustered networks. 
In Fig. 8 we show the dependence of Tco on the triangle 
density v for (fc) = 14 and 18. As expected from the data 
shown in Fig. 7, Tco decreases for rising v, and this de¬ 
pendence turns out to be linear in both cases considered 
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FIG. 7: (a) Crossover temperature as a function of the mean 
degree (fc), for networks with 7 = 3 and A = 8000 nodes. 
Symbols represent results for three v values: 0 (squares), 1 
(circles), and 2 (diamonds). The data shown were obtained 
for networks with several minimum degrees, kg > 3. (b) Ratio 
Tco/J{k) vs (k) for the same kind of networks as in (a). Lines 
are guides to the eye. 


here. The slope of this line is more negative for (fc) = 14 
than for (fc) = 18. In fact, we found: dTco/dv = -0.70 J 
and -0.56 J for (fc) = 14 and 18, respectively. 

This means that for networks with given size N and 
mean degree (fc), including triangles in the networks (i.e., 
increasing the triangle density v) reduces the crossover 
temperature Tco- The reason for this is the following. To 
have a constant (fc) when changing v for a given 7 (7 = 3 
here), one needs a minimum degree kg = {k)/2 — v [see 
Eq. (0}], so that a rise of v is associated to a decrease 
in kg. This causes a reduction of (fc^), and therefore the 
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FIG. 8 : Crossover temperature T^o/J for scale-free networks 
with 7 = 3 and N = 8000 nodes, as a function of the triangle 
density v. Each set of data corresponds to a fixed value of 
the mean degree: (fc) = 14 (circles) and (fc) = 18 (squares). 
Dashed lines are least-square fits to the data points. 


temperature T^o decreases. 

In other words, adding triangles to a network changes 
the degree distribution P(k) itself, apart from introduc¬ 
ing clustering into the network. When one includes trian¬ 
gles {v increases) without changing the other parameters 
defining the networks (i.e., TV, fco, 7), one finds a clear in¬ 
crease in Tco, as shown in Fig. 6 . However, if one changes 
v subject to some particular restriction on the degree dis¬ 
tribution, such as keeping constant the mean value (fc), 
one may find other types of dependence of Tco on v (as 
the decrease shown in Fig. 8 ). 

This suggests that a relevant point here is a compari¬ 
son between clustered and unclustered networks with the 
same degree distribution P{k). This will give insight into 
the ‘direct’ effect of clustering on the critical properties 
of the Ising model. As commented above, the distribu¬ 
tion P{k), as well as Eq. (ITTI) predicting the crossover 
temperature, do not include any information on the clus¬ 
tering present in the considered networks, but only on 
the degrees (connectivity) of the nodes. Thus, a natural 
question is the relevance of the difference between the 
crossover temperature corresponding to networks with 
the same size N and degree distribution P{k), but in¬ 
cluding triangles or not. We have seen above that in the 
case 7 > 3 both types of networks have the same tran¬ 
sition temperature Tc, which agrees with Eq. (fTT|l . This 
is not clear, however, for 7 = 3. To clarify this point 
we have generated networks with the same P{k) as those 
studied above for different v values, but without includ¬ 
ing triangles. In this case, we used the distribution P{k) 
to define the set of degrees employed to build up 

the networks (ki = Si for all i). 

In Fig. 9 we present the temperature Tco as a func¬ 


FIG. 9: Grossover temperature as a function of system size 
N for networks with 7 = 3, iz = 2, and minimum degree ko = 
3. Circles and squares represent results of MC simulation 
for networks with and without triangles, but with the same 
degree distribution P{k). Error bars are less than the symbol 
size. The solid line is the analytical prediction, Eq. (HU, 
obtained from the mean values (fc) and (fc^). 


tion of system size for networks with (circles) and with¬ 
out (squares) triangles, as derived from our MC simula¬ 
tions. Clustered networks were generated with 7 = 3, 
V = 2, minimum degree fco = 3, and several sizes up 
to = 2 X 10® nodes. Unclustered networks were 
built up with the same distribution P{k) as the clus¬ 
tered ones. The results presented in Fig. 9, indicate 
first that clustering increases Tco, and this increase be¬ 
comes more important for larger system size. As dis¬ 
cussed above, for the clustered networks we find a de¬ 
pendence Tco/J = AlnN + B, with a slope A close to 
unity {A = 1.09). We have also plotted the crossover 
temperature predicted by Eq. m, which takes only into 
account the average values (fc) and (/c^) (solid line). A 
linear dependence on In N is also found in this case, but 
with a smaller coefficient A = 0.87. In fact, this line 
crosses with that derived from MC simulations for a sys¬ 
tem size N Ki 10®. 

Looking at the size dependence of the crossover tem¬ 
perature for unclustered networks shown in Fig. 9 
(dashed-dotted line), we observe that for N < 10'^ the 
curve Tco{N) found for these networks is parallel (slightly 
below) to that found for clustered networks (dashed line). 
For larger system sizes, the dashed-dotted curve becomes 
parallel to that corresponding to the analytical model 
(solid line). 

For a given system size, clustered networks display a 
crossover temperature Tco larger than unclustered net¬ 
works with the same degree distribution P{k). This 
means that clustering (triangles in this case) favors an 
increase in Tco, i-e. the FM phase is stable in a broader 
temperature range. This has been already observed in 
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For 2 < 7 < 3, one has for the mean degree: 

{k) « fco ^ , (23) 

7-2 

as in Eq. For 7 = 2, (k) diverges to infinity in the 
large-size limit as 


(k) « - fco In iV . (24) 

Thus, one expects a size-dependent crossover tempera¬ 
ture 


~ ~ 1 _ 1 jy(3-7)/2 

J ^ {k) ^ 3-7 ° 

for 2 < 7 < 3, and 


(25) 


FIG. 10: Crossover temperature T^o/J for scale-free net¬ 
works with i/ = 0 and ko = 3, as a function of system size 
in a logarithmic plot. Symbols represent results for three 
values of the parameter 7 : 2 (squares), 2.5 (circles), and 3 
(diamonds). Error bars are less than the symbol size. Lines 
are guides to the eye. 


the results shown in Fig. 6 , but in that case the difference 
between clustered and unclustered networks was larger, 
due to the inclusion of triangles for > 0 , which changed 
the actual degree distribution P{k) with respect to the 
case 1 / = 0. 

The result for unclustered networks (no triangles) 
shown in Fig. 9 converges to the analytical data given 
by Tco ^ AJ\nN, with A = 0.87, as expected from the 
asymptotic limit for Eq. (ED. Eor large system size, the 
behavior of unclustered networks is in this respect con¬ 
trolled by nodes with large degree. Given that the effec¬ 
tive degree cutoff scales as fccut [see Eq. (IA7I) ]. 

nodes with large k appear progressively as N is increased. 
Thus, the presence of nodes with high degree favors fer¬ 
romagnetic correlations in clustered networks, and there¬ 
fore an increase in the crossover temperature Tco- 


C. Case 7 < 3 

As commented above, for 7 < 3 it is known that corre¬ 
lations between degrees of adjacent nodes appear in scale- 
free networks when no multiple and self-connections are 
allowed, unless one takes a degree cutoff fccut ^ VN [sJ. 
Eor this reason, we generated clustered and unclustered 
networks with 7 < 3 assuming a cutoff kcnt = VN- This 
means that in this case Eq. dS]) does not apply. With 
a calculation similar to that presented in the Appendix, 
one finds in this case for unclustered scale-free networks 
(^ = 0 ): 


(fc2) ~ 1 — 1 ]y(3-A/2 ^ ( 22 ) 

3-7 


Tco ,, ^ VN 

J InN 


(26) 


for 7 = 2 . 

We first present results for unclustered networks {v = 
0). In Fig. 10 we show the temperature Tco as a function 
of system size N for three values of 7 in a logarithmic 
plot, as derived from our MC simulations for networks 
with minimum degree /cq = 3. The exponent 7 increases 
from top to bottom: 7 = 2, 2.5, and 3. For a given system 
size, Tco decreases as 7 rises. This is a consequence of 
a decrease in the ratio (fc^)/(fc) for increasing 7 , which 
causes a reduction in Tco, as predicted by Eq. (ED, where 
the term IvG- 7)/2 dominates for large N. 

For networks with 7 < 3 and large-enough size, log Tco 
derived from the MC simulations is found to display a 
linear dependence on logiV, as expected for a crossover 
temperature diverging as a power of the system size Tco 
with an exponent 2 ; dependent on the parameter 7 
[sslj . Such a linear dependence is obtained for system 
sizes N > No, the size Nq increasing with the exponent 
7 and eventually diverging for 7 —^ 3 , for which Tco 
InA^ is expected (see Fig. 6 ). According to Eq. (ED, 
one expects z = {3 — 'y)l2. In fact, for 7 = 2.5 we find 
2 ; = 0.27 from a fit to the data points corresponding to 
networks with N > 10^. This value of z could further 
decrease for larger system sizes, what is compatible with 
the exponent z = 0.25 expected from Eq. (ED- 

We note that the degree cutoff is relevant for the size- 
dependence of the crossover temperature. Thus, the cut¬ 
off employed here for networks with 7 < 3 (fccut = '/N) 
yields an exponent z = (3 — 7 )/ 2 , to be compared with 
z = (3 — 7 ) 7(7 — 1), given by the “natural” cutoff in 
Eqs. and (IA7|) [see Refs. The latter cutoff 

is known to introduce undesired correlations in networks 
such as those considered here with 7 < 3, as commented 
above. 

We now turn to clustered networks. In Eig. 11 we show 
the crossover temperature vs system size N for networks 
with different triangle densities: = 0, 1, and 2. In 

panels (a) and (b) results are given for 7 = 2 and 2.5, 
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FIG. 11: Crossover temperature Tco/J as a function of sys¬ 
tem size N for scale-free networks with minimum degree ko 
= 3, and three values of the triangle density v: 0 (squares), 
1 (circles), and 2 (diamonds), (a) Networks with 7 = 2; (b) 
networks with 7 = 2.5. Lines are guides to the eye. 


respectively. For small network size, Tco appears to be 
higher for larger parameter v in both cases. This differ¬ 
ence, however, decreases as N is increased, and for each 
value of 7 the results for different triangle densities iz 
converge one to the other. Thus, in the logarithmic plots 
of Fig. 11 differences between the crossover temperature 
for different v values become irrelevant for system size 
N - 10®. 

Similarly to the case 7 = 8 presented in Sec. IV.B, for 
7 < 3 it is also interesting to compare results for clustered 
and unclustered networks with the same degree distribu¬ 
tion P{k). In Fig. 12 we present the temperature Tco as 
a function of system size for networks with (circles) and 



FIG. 12: Crossover temperature as a function of system size 
N for networks with 7 = 2, iz = 2, and minimum degree ko = 
3. Circles and squares represent results of MC simulations for 
networks with and without triangles, respectively, but with 
the same degree distribution P{k). Error bars are less than 
the symbol size. The solid line indicates the result derived 
from {k) and (A:^) by using Eq. dill) . 


without (squares) triangles, as derived from our MC sim¬ 
ulations. Clustered networks were generated with 7 = 2 , 
V = 2, minimum degree ko = 3, and various system sizes. 
For N < 2000 nodes, results for clustered and unclus¬ 
tered networks coincide one with the other within statis¬ 
tical noise. For N > 2000, both sets of data progressively 
separate one from the other, so that the crossover tem¬ 
perature for clustered networks (including triangles) is 
larger than that corresponding to unclustered networks 
(no triangles). The solid line in Fig. 12 is the analyti¬ 
cal prediction obtained from Eq. m by introducing the 
mean values (k) and {k'^) corresponding to the actual 
networks. We observe something similar to the case of 
7 = 3 shown in Fig. 9. For large N, results of MC sim¬ 
ulations for unclustered networks approach the analyti¬ 
cal expectancy, lying below the solid line, whereas data 
for clustered networks become higher than the analytical 
prediction, and progressively deviate from the latter as 
system size N is increased. 

As for the case 7 = 3, we conclude that clustering fa¬ 
vors a stabilization of the FM phase vs the PM one, hence 
increasing the crossover temperature Tco- This becomes 
more noticeable for larger network size, where nodes with 
higher degree progressively appear. 


V. SUMMARY 

We have studied the FM-PM transition for the Ising 
model in clustered scale-free networks by means of Monte 
Carlo simulations, and the results were compared with 
those found for unclustered networks. Our results can be 
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classified into two different regions, as a function of the 
exponent 7 defining the power-law for the degree distri¬ 
bution in scale-free networks. 

For 7 > 3, we find in all cases a well-defined tran¬ 
sition temperature Tc in the thermodynamic limit, in 
agreement with earlier analytical calculations and MC 
simulations. This refers equally to clustered and unclus¬ 
tered networks. Adding motifs (triangles in our case) to 
the networks causes an increase in the transition tem¬ 
perature, as a consequence of the associated change in 
the degree distribution P(/c), and in particular in the 
mean values {k) and (A:^). However, for clustered and 
unclustered networks with the same degree distribution 
P{k), one finds no difference in Tc, which coincides with 
that predicted by analytical calculations [ Eq. (fTTll ] . This 
conclusion agrees with that drawn in Ref. l5lL where the 
addition of motifs was found to increase the transition 
temperature, without changing the critical behavior. 

For networks with 7 < 3, the situation is different. In 
this case, the crossover temperature Tco increases with 
system size N. For 7 = 3 we found Tco JlnA^, and 
for 7 < 3 we obtained Tco ^ JN^, with an exponent 
z = (3 — 7 )/ 2 . Comparing clustered and unclustered 
networks, the conclusions obtained for 7 < 3 differ from 
those found for 7 > 3. 

For 7 < 3, Tco is similar for clustered and unclustered 
networks with the same degree distribution P(k), when 
one considers small network sizes {N < 10^). This behav¬ 
ior changes for larger networks, and the crossover tem¬ 
perature Tco of clustered networks becomes progressively 
larger than that corresponding to the unclustered ones. 
Thus, we hnd that FM correlations are favored by includ¬ 
ing triangles in the networks, in particular in the presence 
of nodes with large degree k. 

It is important to note that this conclusion refers to a 
comparison of clustered and unclustered networks with 
the same degree distribution P{k). Caution should be 
taken when comparing results for networks with different 
degree distributions. Thus, when triangles are added on 
a network with a given distribution of links, the crossover 
temperature Tco increases [see Fig. 6 ], but in this case the 
rise in Tco is mainly due to the change in P{k) caused 
by the inclusion of triangles. However, for clustered net¬ 
works with the same exponent 7 (large-degree tail) and 
mean degree (fc), an increase in triangle density v causes 
a decrease in Tco [see Fig. 8 ], as a consequence of the 
associated variation in the distribution P{k), and in par¬ 
ticular in the mean value (A;^). 

We also note that results for unclustered or clustered 
scale-free networks may differ when different degree cut¬ 
offs are employed, especially for 7 < 3, as the depen¬ 
dence of their properties on system size can be appre¬ 
ciable. This applies, in particular, to the scaling of the 
crossover temperature Tco on N for large system size. To 
avoid correlations between degrees of adjacent nodes, we 
have employed here for power-law distributions a degree 
cutoff A:cut = However, for clustered networks the 

presence of triangles introduces degree correlations. 


Other distributions different from the short-tailed 
Poisson-type introduced here for the triangles, could be 
considered to change more dramatically the long-degree 
tail of the overall degree distribution P{k). For exam¬ 
ple, a power-law distribution for the triangles (with an 
exponent 7 ') may give rise to an interesting competition 
between the exponents of both distributions (for single 
links and triangles), which could change the critical be¬ 
havior of the Ising model on such networks as compared 
to those discussed here. 
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Appendix A: Finite-size effects on degree 
distribntions 


Here we present some expressions related to finite-size 
effects in scale-free networks. We assume a degree distri¬ 
bution Psf{k) defined as 



nk 

0 , 


for k > ko 
for k < ko 


(Al) 


with 7 > 2 and n a normalization constant. 

Due to the finite network size N, an effective cutoff 
fcfii t. app ears for the degree distribution in these networks 
[ 3 ^ l3^ . This cutoff is such that ^ Psf(fc) ^ l/A^j 
indicating that the number of nodes with k > kcut is 
expected to be on the order of unity. For concreteness, 
we write 


J2PAk) = ^ (A2) 

^cut 


with c = 0{1). Replacing the sum by an integral, we find 


■ (A3) 

7 — i 

^cut 


For c/N <C 1, one has for the normalization condition: 


^cut 




(U-’ - . (A4) 


kg 


7-1 


Then, the normalization constant n is: 


n = 


~ (7 - 1) ^0 ^ ■ 


jLi -7 _ l.k-1 
'*'0 ''"cut 


Combining Eqs. 


, and 
7-1 


one finds 


* A:o * 


(A5) 


(A 6 ) 
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and for fccut ^ ko: 


3, we find: 


N\ "'-1 


ko 


so that fccut as in Refs, 

this cutoff, one has for the mean degree 


(A7) 

Considering 


fccut -1 7 2—"7 7 2—7 

(fc) = ^ A: (fc) = ^ 


fco 


’ '^cut '^0 


7 _ 7 „ 1-7 


(A8) 


which gives, using Eq. (TATI) : 


(fc) Ri (fc)c 


-V — ^ 


(A9) 


with (fc)oo ~ fco (7 — l )/(7 ~ 2). Thus, we find for the 
mean degree (k): 


{k^) « {k^)c 




(v)' 


with (fc^)oo ~kl{'^ - 1 )/(7 - 3). 
For 7 = 3, we have 


2\ 


= n 


1 


In ■ 


^ k' fco 2 _ fc ^^2 - 


and using Eq. (I ATE we obtain for fccut ^ k^: 

(fc 2 ) = kl In 7+ 0 ( 1 ) 


(All) 


(A12) 


(A13) 


(fc) 


{k} 00 




(AlO) 


A similar calculation can be carried out to estimate 
finite-size effects on (fc^) for scale-free networks. For 7 > 


For 7 < 3 a strict cutoff fccut = '/N has been intro¬ 
duced in the networks discussed in the present paper, 
in order to avoid undesired correlations [^. Thus, the 
equations given in this appendix do not apply to the ac¬ 
tual networks discussed in Sec. IV.C. 
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